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ABSTRACT 

We report first-epoch circular polarization results for 133 active galactic nuclei in the MOJAVE program to 
monitor the structure and polarization of a flux limited sample of extra-galactic radio jets with the VLBA at 15 
GHz. We found strong circular polarization (> 0.3%) in approximately 15% of our sample. The circular polar- 
ization was usually associated with jet cores; however, we did find a few strong jet components to be circularly 
polarized. The levels of circular polarization were typically in the range of 0.3-0.5% of the local Stokes-/. We 
found no strong correlations between fractional circular polarization of jet cores and source type, redshift, EGRET 
detections, linear polarization, or other observed parsec-scale jet properties. There were differences between the 
circular-to-linear polarization ratios of two nearby galaxies versus more distant quasars and BL Lac objects. We 
suggest this is because the more distant sources either have (1) less depolarization of their linear polarization, 
and/or (2) poorer effective linear resolution and therefore their VLBA cores apparently contain a larger amount of 
linearly polarized jet emission. The jet of 3C 84 shows a complex circular polarization structure, similar to obser- 
vations by Homan & Wardle five years earlier; however, much of the circular polarization seems to have moved, 
consistent with a proper motion of 0.06c. The jet of 3C273 also has several circularly polarized components, and 
we find their fractional circular polarization decreases with distance from the core. 

Subject headings: galaxies : active — galaxies: jets — quasars : general — radio continuum: galaxies — BL 
Lacertae objects: general — polarization 



1. INTRODUCTION 

Only a small fraction of the radio emission from extra- 
galactic radio jets associated with Active Galactic Nuclei 
(AGN) is circularly polarized; however, this weak component 
is a potentially important probe of the particle population and 
magnetic field structure of the jets. Circular polarization (CP) 
can be produced either as an intrinsic component of the emitted 
synchrotron radiation, or via Faraday conversion of linear po- 
larization into circular (Jones & O'Dell 1977). If the intrinsic 
mechanism dominates, jets must be a predominantly electron- 
proton plasma and must contain a significant fraction of uni- 
directional ordered magnetic field. If the conversion mecha- 
nism dominates, circular polarization probes the low energy end 
of the relativistic particle distribution (e.g. Wardle & Homan 
2003), a key parameter in studying the bulk kinetic luminosity 
of AGN jets and their particle content (Celotti & Fabian 1993; 
Wardle et al. 1998). Faraday conversion is expected to be the 
dominant mechanism on theoretical grounds (e.g. Jones 1988); 
however, there is still only limited observational evidence sup- 
porting this view (Wardle et al. 1998; Homan & Wardle 2004). 

Early integrated measurements of CP from AGN jets (Weiler 
& de Pater 1983; Komesaroff et al. 1984) revealed circular po- 
larization to be only a very small fraction, typically < 0.1%, 
of their total flux. However, the Very Long Baseline Array 
(VLBA)', has made it possible to study circular polarization 
from extragalactic jets at sub-milliarcsecond resolution. War- 
dle et al. (1998) and Homan & Wardle (1999, hereafter HW99) 
found circular polarization in cores of four powerful AGN jets 

' The VLBA is operated by the National Radio Astronomy Observatory, a facility 
Associated Universities, Inc. 
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at local fractional levels of 0.3 to 1% of Stokes-/ with the 
VLBA. Circular polarization has since been detected with vari- 
ous instruments in other AGN jets (Rayner et al. 2000; Homan 
et al. 2001) and a wide variety of other synchrotron emitting 
sources: low-luminosity AGN such as Sagittarius A* and M81* 
(Bower, Falcke, & Backer 1999; Sault & Macquart 1999; Brun- 
thaler et al. 2001; Bower, Falcke, & Mellon 2002), intra-day 
variable sources (Macquart et al. 2000), and galactic micro- 
quasai-s (Fender et al. 2000, 2002). 

Despite the detection of the circular polarization in a vari- 
ety of sources, observations of circular polarization from AGN 
jets are fragmented and incomplete. A large temporal gap ex- 
ists between the integrated measurements of the 1970's to early 
1980's and the high resolution VLBA measurements from the 
mid 1990s onward. Only a few AGN jets have high resolution, 
multi-epoch circular polarization observations that are matched 
in observing frequency and sensitivity. An even smaller num- 
ber of sources have multi-frequency circular polarization ob- 
servations suitable for measuring a spectrum. Additionally, the 
small, ad hoc samples studied to date have not allowed for an 
analysis of correlations between the appearance of strong cir- 
cular polarization and other source properties. While these pre- 
vious studies raised a number of interesting possibilities of how 
circular polarization could be used to constrain the magnetic 
field properties and particle populations of AGN jets, a system- 
atic survey of a large complete sample is necessary to make 
further progress. 

We have exploited the VLBA's unique imaging capabilities 

of the National Science Foundation operated under cooperative agreement by 
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by studying the milliarcsecond-scale circular polarization of 
a complete AGN jet sample at 15 GHz as part of the MO- 
JAVE (Monitoring Of Jets in AGN with VLBA Experiments) 
program. The MOJAVE sample of 133 AGN is complete and 
flux limited on the basis of VLBA-scale emission. The selec- 
tion criteria are described in detail in Lister & Homan (2005) 
which also reports first-epoch linear polarization results; how- 
ever, in brief, for inclusion in MOJAVE, a source must have 
had a total 15 GHz VLBA flux density of at least 1.5 Jy (> 2 
Jy for sources south of the celestial equator) at any epoch from 
1994-2003. When completed, MOJAVE will have obtained at 
least four epochs of full Stokes VLBA images at 15 GHz of all 
sources in the sample. In addition, multi-frequency follow-up 
programs will study the circular polarization spectrum of the 
strongly circularly polarized sources identified by MOJAVE. 

In this paper, we report only the first-epoch circular polar- 
ization results from the MOJAVE program. Future papers will 
explore other aspects of this rich data set including multi-epoch 
circular polarization variability. The details of our observations 
and calibration appear in section 2. We present and discuss our 
results in section 3, and our conclusions appear in §4. 

2. OBSERVATIONS AND CALIBRATION 

Each MOJAVE observing session targets 18 sources over a 
24 hour period, and this paper presents results drawn from the 
first eleven sessions, from 2002 May 31 through 2004 February 
11. During each 24 hour session, observations of the 18 sources 
are broken up into highly interleaved scans, each approximately 
seven minutes in duration. The schedule is designed to maxi- 
mize coverage in the (m, v)-plane; however, as discussed below, 
the highly interleaved nature of the observations also improves 
our ability to calibrate the antenna gains to detect circular po- 
larization. We observed each source for approximately 65 min- 
utes of total integration time. No additional calibrators were 
observed as all of our target sources were strong and compact. 

The data were recorded at each VLBA antenna at a cen- 
tral observing frequency of 15.366 GHz using 1-bit sampling 
and were correlated at the VLBA correlator in Socorro, NM. 
The correlator output was 2 s integrations for all four cross- 
correlations (RR, RL, LR, LL). The correlated data contained 
four intermediate frequencies (IFs) and 16 channels per IF, 
yielding a total bandpass of 32 MHz in both RR and LL. 
Tapes of the correlated data were shipped to NRAO in Char- 
lottesville, VA and Purdue University where they were loaded 
into NRAO's Astronomical Image Processing System (AIPS) 
for further processing. 

The basic calibration procedure followed the standard meth- 
ods outlined in the AIPS cookbook and the procedure by HW99 
for measuring circular polarization with the VLBA using the 
gain transfer technique. Below we briefly describe our basic 
calibration procedure, in §2. 1 we discuss some improvements 
to the gain transfer procedure developed for these observation, 
and in §2.2 we discuss some checks on our calibration proce- 
dure. 

The data were loaded into AIPS using the VLBALOAD pro- 
cedure. The procedure VLBACALA was then used to do a- 
priori amplitude calibration, running both the ACCOR and AP- 
CAL tasks, including an atmospheric opacity correction. The 
parallactic angles were removed using the procedure VLBA- 
PANG and the pulse calibration was applied using VLBAP- 

^ For the total intensity and linear polarization imaging described in Lister & I 

= Imodel and LL = Imodel- 



COR. The procedure VLBAFRNG was then used to perform 
fringe fitting. All of our sources were strong enough to serve 
as their own fringe calibrator, and no special interpolation from 
strong to weak sources was necessary. We used the procedure 
CRSFRING to remove any residual multi-band delay difference 
between the right and left-hand systems of the array. Finally 
we ran BRASS to apply a bandpass correction before averaging 
across the channels within each IF and 'splitting' the data off 
into individual source files. 

Following a 2 s time-scale point-source phase calibration to 
increase coherence, the data were written out to the CALTECH 
VLBI program, DIFMAP (Shepherd 1997). In DIFMAP, tiie 
data were averaged to 30 seconds and edited in a station-based 
manner. An edited version of the data was then saved with- 
out any self-calibration applied. We will refer to this ver- 
sion of the data as the edited/unself-calibrated data. Normal 
self-calibration and imaging techniques were then used within 
DIFMAP to obtain the best Stokes-/ model for each source. 

For each source, the best Stokes-/ model and the 
edited/unself-calibrated data were then read back into AIPS, 
and the data were self-caUbrated against the model using the 
task CALIB in amplitude and phase (SOLMODE = 'A&P') 
on a 30 second timescale. The antenna feed-leakage term cor- 
rections from LPCAL (see below) were applied using SPLIT 
and the parameter DOPOL = 3 to guarantee that the correc- 
tions were applied to both the parallel and cross-hand corre- 
lations. A final pass of CALIB self-calibrated the data again 
against the Stokes-/ model to remove any errors by the initial 
self-calibration prior to D-term removal. For circular polariza- 
tion calibration via gain transfer it is important that both passes 
of CALIB described above were performed only making the 
rigorous assumption^ that {RR + LL)/2 = Imodel- This left as un- 
calibrated any gain offsets between the right (R) and left (L) 
hand feeds at each antenna. The gain transfer calibration de- 
scribed in §2.1 was used to remove these offsets. After correc- 
tion of the R/L antenna gains, final images were made in all 
four Stokes parameters. 

It is important to note that the antenna feed leakage correc- 
tions were initially determined on each source individually us- 
ing LPCAL, and the final corrections (which were applied to 
the data) were determined by taking the median of these af- 
ter removing obvious outUers. Our typical feed leakage terms 
(D-terms) had amplitude of about 0.02 with the largest being 
0.05, and we estimate that we have made these corrections 
accurate to about ±0.005 . Following HW99, we estimate 
any uncertainty in our fractional circular polarization measure- 
ments, nic, due to leakage correction errors to be of the order 
Amc < O.OOIwl-i- 1 X 10"^, where ml is the fractional linear 
polarization of the map peak (typically a few percent or less). 
Thus any uncertainty due to leakage feed corrections is very 
much smaller than the uncertainty due to the gain calibration 
and can be neglected. 

2.1. R/L Gain Calibration 

As described in HW99 and Homan et al. (2001), the gain 
transfer technique requires that we initially make no assump- 
tion about circular polarization during self-calibration, i.e., we 
only assume {RR+LL)/2 = Imodel as discussed above. This will 
leave as un-calibrated a small relative R/L complex gain ratio at 
each VLBA anterma. The essential idea is to calibrate this gain 

(2005) we assumed there was no Stokes-V in these CALIB steps by assuming 
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ratio with a smoothed set of R/L complex antenna gains deter- 
mined from all 18 sources in each 24 hour observing session. 

To determine the smoothed set of antenna gains that will ulti- 
mately be applied to the data, we first solve for the R/L complex 
gain ratios by self-calibrating each source assuming no circular 
polarization in that source (RR = Imodei and LL = Imodei)- These 
solutions are not applied to the data, but rather they are col- 
lected into a single caUbration table which we smooth by using 
a six-hour running median with all sources equally weighted. 
Solutions from sources with low gain SNR or very strong cir- 
cular polarization (greater than 0.5%) are prevented from con- 
tributing to the running median^. To account for the possibil- 
ity of smaller levels of real circular polarization showing up in 
the raw R/L antenna gains of individual sources, we subtract 
from the raw R/L gain ratios for each source a single mean off- 
set (taken across aU times, antennas, and IPs for that source) 
from the running median. The sum of the applied offsets for 
all sources in a given 24 hour observing session is constrained 
to be zero. This constraint guarantees quick convergence of 
an iterative algorithm which generates a trial set of smoothed 
gains, computes mean offsets for each source, takes those off- 
sets out of the data, and then repeats the process. The constraint 
also underscores the fundamental assumption of the gain trans- 
fer technique: that any single observing session will have an 
average circular polarization of zero across aU the sources in 
the sample"*. This assumption may introduce a small overall 
gain bias in any individual epoch; however, our Monte Carlo 
method of estimating the uncertainty due to the gain calibration 
(described below) naturally accounts for this possibihty. 

We note that the procedure described above for accounting 
for sources with small levels of real circular polarization differs 
from the procedure used in HW99 and Homan et al. (2001). 
In those papers, we prevented such sources from contributing 
to the smoothed gains through an iterative process where the 
strongest polarized sources were removed first, the smoothed 
gains recomputed, and the process repeated until no source con- 
tributing to the smooth gains appeared to have real circular po- 
larization. In those papers, we found the final results were in- 
sensitive to the detailed steps of this process; however, the pro- 
cess was time consuming and difficult to automate. For the 
MOJAVE program, we needed a procedure that was automated, 
reproducible, and could be easily tested by a Monte Carlo sim- 
ulation we developed. As a check that the new procedure gave 
the same results as our previous methods, we reprocessed the 
antenna gain solutions from the 40 source observation analyzed 
in Homan et al. (2001). For the circularly polarized sources 
in that experiment, we found an RMS R/L gain difference of 
0.017% between our technique and the original. This difference 
is approximately one-third of the estimated gain uncertainty of 
that experiment, verifying that the two techniques indeed give 
effectively the same results. 

2.1.1. Uncertainty in Measured CP due to Gain Transfer 

In the appendix of Homan et al. (2001), we laid out a pre- 
scription for estimating the uncertainty in circular polarization 
measurements made by gain transfer Here we improve on that 
work by using a Monte Carlo simulation to test our methods 
and obtain more accurate uncertainty estimates. As in Homan 
et al. (2001), we focus on amplitude gain errors, as those domi- 

^ Such strongly polarized sources were very rare; however, when we did find one, 
contributing to the median gain solution 

AGN as a whole have no preferred sign of circular polarization, and thus we expect 



nate the uncertainty in the technique. Phase errors will not con- 
tribute to first order for point-like sources, and for sources with 
significant extended structure, phase errors will produce anti- 
symmetric structure around the phase center. In section §2.2 
we discuss possible phase effects and tests to constrain/resolve 
those issues. 

The Monte Carlo simulation we developed generates 18 
point sources (Most sources in our sample are dominated by 
a point source, e.g.. Lister & Homan (2005).), and for each 
source generates a circular polarization with magnitude drawn 
at random from the distribution of observed mean gain offsets 
in the real gain table. The sign (-t- or -) of the CP signal for 
each source is determined at random. Finally, a Gaussian devi- 
ate representative of our measurement uncertainty is added to 
each value. Each of these generated CP values is then stored 
for comparison to the "measured" value at the end of a Monte 
Carlo trial. 

To generate the gain table for the simulation, a random 
slowly-varying gain curve is created for each antenna/lF with 
amplitudes and timescales similar to those in the real gain 
curves. To this slowly-varying gain curve is added randomly- 
generated scan-to-scan variations of the same range of ampli- 
tudes as seen for the corresponding source in the real gain table. 
Finally, the generated CP values are added to the gain curves as 
offsets for each source. This step is important to mimic the 
real gain table, which was generated from the assumption that 
each source has no circular polarization, so the real raw gain 
table contains 'corrections' to the circular polarization of each 
source in addition to the actual gain variations. 

Finally, this randomly-generated 'raw' gain table is pro- 
cessed by exactly the same procedure we used on the real gain 
curve. A smoothed version of the raw table is produced by our 
standard procedures and compared with the raw gain table. If 
the procedures work correctly, the only differences between the 
raw and smoothed tables should be the circular polarizations set 
at the beginning of the simulation trial. Therefore, by compar- 
ing the two tables, we can estimate the uncertainty due to the 
gain calibration procedure. We repeat this simulation for 1000 
independent trials to generate an estimate of the uncertainty for 
each source. 

To safe-guard against the possibility that an unusual combi- 
nation of circularly polarized sources in our data has made our 
results more imcertain than the above simulation indicates, we 
also ran a second version of the Monte Carlo described above. 
The only difference is that the initial circular polarization of 
each source in the simulation is forced to match the mean gain 
offset of the corresponding source in the real gain table (plus 
a Gaussian deviate to represent measurement uncertainty). To 
these values, we add a randomly-determined overall bias (the 
same for all 18 sources) to account for the possibility that the 
circular polarizations in the sample do not average to zero. The 
size of this gain bias is determined from the RMS of the mean 
offset values in the real gain table. The larger of the uncer- 
tainties from the first and second version of the Monte Carlo 
simulations are used to represent the total gain uncertainty for 
a given source; however, in practice the two Monte Carlo sim- 
ulations give very nearly the same estimate of uncertainty. 

Figure 1 shows a plot of uncertainty in fractional circular po- 
larization measurements due to gain calibration as a function 

we simply repeated the caUbration with the source in question prevented from 
the net circular polarization from a large sample of AGN to be close to zero. 
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of the peak Stokes-/ flux for the 133 sources in the MOJAVE 
sample. The trend of larger uncertainties for weaker sources 
is expected and is driven by the larger self-caUbration gain un- 
certainties in the weaker sources; however, it will be important 
to consider this trend when interpreting the results in §3. Note 
that most sources in the sample have uncertainty contributions 
from gain calibration very close to 0. 1 % , and given that the gain 
uncertainties dominate the total uncertainty (the much smaller 
RMS map noise is added in quadrature), our procedures are able 
to robustly detect CP signals as small as 0.3% of the Stokes-/ 
peak at the 3-a level. 

2.2. Phase Calibration for Circular Polarization 

Because there is no phase connection between sources in our 
sample (they are well separated in the sky and fringe-fit sepa- 
rately), the gain caUbration technique described above is only 
effective for caUbrating the anterma amplitude gains. We do 
allow the possibility of a phase correction in the gain transfer 
technique (the median phase difference between the right and 
left hands is treated in a similar fashion to the median amphtude 
ratio.); however the applied corrections are generally < 0.1 de- 
grees for each antenna/IF combination. 

For effective relative R/L phase correction at each antenna 
we rely first on the pulse calibration corrections, then on fringe 
fitting and the initial phase cahbration in CALIB to increase co- 
herence before averaging the data in time in DIFMAP. Each of 
these steps provides separate phase corrections for the left and 
right-hand feeds at each antenna. The pulse calibration makes 
no assumptions about circular polarization in the source; how- 
ever both fringe fitting and the initial phase caUbration to in- 
crease coherence make their corrections assuming the source is 
point-like. Any real circular polarization on a point source will 
not appear as a gain error in the antenna phases (HW99), so 
for most sources in our sample phase corrections found in this 
fashion will remove only genuine R/L phase gain errors from 
the data. 

If a source has significant extended structure in Stokes-/, the 
right and left hand systems of the array may receive shghtly dif- 
ferent solutions in the algorithm's attempt to find point source 
solutions to the right and left hand data separately. Most likely 
this will generate a small difference in the phase centers of the 
right and left hand data. While this small difference will not 
affect the Stokes-/ image, it may be enough to generate an anti- 
symmetric stripe of circular polarization across the strongest 
parts of the source. The effect is only apparent in the circu- 
lar polarization image because Stokes-V is formed by taking 
the difference of flie right and left hand data {V = (RR-LL) /2), 
whereas Stokes-/ is formed from their average. See panels (a) 
and (c) of Figure 2 for examples of this effect. 

To account for this effect and other possible phase errors in 
our data, we have generated two circular polarization images of 
each source. The first image is generated from the gain transfer 
calibration alone and may contain unaccounted-for phase er- 
rors. The second image is generated after an additional round of 
self-calibration of the phases alone. During this self-calibration 
of the phases, we assume there is no circular polarization in 
the data so both the RR and LL data are compared to Imodei to 
find the antenna phase corrections for the right and left-hand 
feeds separately. This self-calibration will very effectively re- 
move any genuine R/L phase gain errors remaining in the data, 
as shown in panels (b) and (d) of Figure 2. The large major- 
ity of the sources in our sample are core-dominated and thus 



show no significant differences between the images produced 
before and after this phase self-calibration; however, an impor- 
tant minority of sources with significant extended structure do 
show differences (usually both in their noise levels and in their 
apparent circular polarization) between the two sets of images. 

Figure 3 illustrates two cases where the antenna phase er- 
rors have been removed to reveal real circular polarization. It is 
important to know to what extent (if any) the resulting circular 
polarization in these images has been modified by the phase cal- 
ibration assuming zero circular polarization. In general, phase 
differences between RR and LL due to real circular polariza- 
tion do not show up in a consistent antenna-based fashion (as 
would phase gain errors); however it is possible that some com- 
bination of antenna-based phase gains could be found by the 
self-calibration algorithm which minimizes these real phase dif- 
ferences and thus modifies the real circular polarization of the 
source. 

We have tested for this effect by (1) removing all circular 
polarization from a real source (by setting RR = LL = I), (2) 
adding known amounts of circular polarization, and then (3) 
self-calibrating the phases of that data assuming no circular po- 
larization. We have done this for several sources from our ob- 
servations, and in general we find that the known circular po- 
larization is not modified very much by this procedure. We 
usually see some reduction in amplitude from a few percent up 
to about 10% of the known CP value, and the location of the CP 
peak may be shifted a fraction of a beam-width to better align 
with the map peak (assuming it is near the map peak to begin 
with). In some sources with a broad core-region (core + barely 
resolved component), a single component of CP is spread out 
by this self-calibration to encompass the whole region. Large 
shifts in CP location (more than a half a beam-width) or large 
reductions in amplitude (more than 20%) do not seem possi- 
ble with this procedure. For example, our tests on the images 
in Figure 3 show that real circular polarization that began on 
the first jet component would remain there after phase self- 
calibration and would not be shifted to the core; thus, we can 
be confident that the real circular polarization in this epoch of 
3C 273 is actuaUy associated with the core as revealed in panel 
(b). 

3. RESULTS AND DISCUSSION 

Our first-epoch circular polarization measurements are sum- 
marized in Tables 1, 2, and 3. Table 1 Usts the linear and cir- 
cular polarization measurements at the peak of each map. In 
most cases the Stokes-1 map peak corresponds to the unresolved 
VLBI core; however, in a handful of cases the core and map 
peak do not coincide. In most of these cases, the peak appears 
at the location of a strong component within the jet, and we 
take the core to be the compact component at the base of the 
jet. The core of 0923+392 is taken at the location identified by 
Alberdi et al. (1993). The core of 0316+413 (3C 84) is a special 
case and is discussed further in §3.3. The off-peak jet cores are 
listed in Table 2, along with measurements of their linear and 
circular polarization at the core location. Table 3 lists the cir- 
cular polarization measured on all jet components in our sam- 
ple with peak fluxes stronger than 300 mJy beam"'. We chose 
300 mJy beam"' as a reasonable cutoff for the jet component 
analysis as our typical uncertainty levels would still allow very 
strong circular polarization of Wc = 0-5% to be detected on such 
a component. 

The linear polarization was always measured at the location 
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of the Stokes-I peak (or at the location of the VLBI core or jet 
component in the cases of Tables 2 and 3 respectively). Limits 
are given for the Unear polarization when the measured value 
does not exceed five times the map RMS noise. Circular polar- 
ization was also measured at the location of the Stokes-1 peak 
(or at the location of the VLBI core or jet component for Ta- 
bles 2 or 3) unless there was a CP peak within two pixels of 
the Stokes-I measurement location. In that case, the value of 
the CP peak was used instead. This procedure was used be- 
cause with our phase calibration (see section 2), there may be 
a small positional shift in the location of the circular polariza- 
tion relative to the Stokes-I. The raw circular polarization mea- 
surements and their uncertainties are given in Tables 1, 2, and 
3. The uncertainties includes both the gain transfer uncertainty 
from our Monte Carlo simulation (see §2) and the RMS noise^ 
added in quadrature. It is important to note that the fractional 
circular polarization values reported in the Tables include both 
measurements and limits. While we do actually measure a CP 
value in every case, most of these measurements are not signif- 
icant, and thus we have chosen a 2a cutoff, below which we 
report a Umit on the value of |mc | . The quoted limits on |mc | are 
the larger of (| V| + Ictv) or 2(Tv. The raw values of rUc in the 
cases that are less significant than 2(Tv can be computed from 
the y and / intensities reported in the Tables. 

We note that we have applied exactly the same procedures to 
the cores and jet components in Tables 2 and 3 as we applied 
to the map peak values; and, as argued in §2.1, the interpreta- 
tion of the gain uncertainties is most robust for the map peak. 
Therefore, we are somewhat less confident of the uncertainties 
reported in Tables 2 and 3; however, we believe that the gain 
uncertainties will still accurately account for the possibility of 
any overall gain bias in the CP map. The gain uncertainties 
do not account for the possibihty of random (but not net) gain 
errors creating small amounts of spurious circular polarization 
off the map peak; however, the RMS noise level in the map 
naturally takes this into account. We notice no particular trend 
for map noise to be scattered preferentially onto weaker source 
structure, and we therefore believe that the combination of gain 
uncertainties and map RMS noise provides a reasonable uncer- 
tainty estimate for the off-peak circular polarization reported in 
Tables 2 and 3. 

3.1. Circular Polarization on VLBI Cores 

Figure 4 shows the distribution of core circular polarization 
for all the sources in our sample. For the purposes of this figure 
we use measured values rather than limits in every case. We 
detect circular polarization in 17 out of 133 jet cores at the 3a 
level or higher significance. Of the 115 sources with overall 
uncertainties less than 0.15%, we measure circular polarization 
stronger than 0.3% of the local Stokes-I in 17 of those sources 
(all but two of those at > 3(7 significance). The existence of 
strong circular polarization (> 0.3%) in about 15% of our sam- 
ple is somewhat higher than the rate of strong circular polariza- 
tion observed by Homan et al. (2001) with the VLBA at 5 GHz. 
They found CP stronger than 0.3% in only 2 out of 36 sources 
(6%); however, these rates of strong circular polarization detec- 
tion are compatible according to Poisson statistics. 

We examined the circular polarization distribution between 
source types (using the optical classifications given by Lister & 
Homan 2005), and a Gehan's GeneraUzed WUcoxon test from 

' The RMS noise was calculated over a large number of pixels away for source structure. For each pixel, we used largest value of CP within 2 pixels of that location. 
This procedure increases our estimate of the RMS noise, but it is consistent with the procediffe we use for measimng CP from a Stokes-I peak. 



the ASURV package (Lavalley et al. 1992) shows no signifi- 
cant difference between the distributions for quasars and BL 
Lac objects. Unfortunately there are simply too few galaxies in 
the MOJAVE sample for a robust test of their distribution. A 
Gehan's test also shows no significance difference between the 
EGRET and non-EGRET detected sources (as defined in Lister 
& Homan 2005) in terms of their CP levels. 

Due to the large number of limits in our sample, we used a 
Kendall-Tau test for censored data (Lavalley et al. 1992) to test 
for correlations between circular polarization and various other 
source properties. Table 4 collects the results of these tests. The 
Kendall-Tau test gives the probability that no correlation exists 
between fractional circular polarization, \mc\, and the property 
in question. Small values in the Table indicate that a correlation 
may exist, and we take p < 0.05 as suggestive of a correlation. 
Of the twenty properties tested in Table 4, we find only three 
with p < 0.05: 15 GHz jet flux (5^^,), 15 GHz core/jet flux ra- 
tio (Rjet), and 15 GHz core/total flux ratio (Rtoiai) as defined 
by Lister & Homan (2005). These three properties are closely 
related through the prominence of the VLBI scale jet, so we 
really only have one property with a possible correlation at the 
p < 0.05 level, consistent with pure chance. Removing 3C 273 
from consideration raises p to 0.12 for Sj^i and 0.09 for R,oiai, 
so we don't find strong evidence that circular polarization is 
correlated with any of these properties. 

While Galactic latitude, \b\, is not correlated with circular 
polarization, it is interesting that of the 12 sources near the 
galactic plane (\b\ < 10 degrees), none of them have circular 
polarization at > 2cr. Treating these sources and those with 
\b\ > 10 degrees as separate samples, a Kolomogorov-Smimov 
gives only a 5.3% chance they are drawn from the same dis- 
tribution. Again, this low probability may be spurious as there 
seems little reason the galactic plane should interfere with cir- 
cular polarization signals. In particular, any galactic Faraday 
rotation should have no noticeable effect on circular polariza- 
tion (e.g., Homan et al. 2001). 

One property that might naturally be expected to correlate 
with core circular polarization is core linear polarization. In 
the most common models of circular polarization production 
(Wardle & Homan 2003), both linear and circular polarization 
increase with increasing field order; however, the two may not 
be correlated in a few scenarios: (1) if circular polarization is 
produced through a high internal Faraday depth model (e.g., 
Ruszkowski & Begelman 2002; Beckert & Falcke 2002), (2) 
if the linear polarization is externally depolarized (e.g., Homan 
et al. 2001), or (3) if the circular and linear polarization we 
observe are produced in two distinct regions of the jet but are 
unresolved in what we see as the VLBI core. 

Figure 5 is a plot of absolute fractional core circular polariza- 
tion, \mc\, versus fractional core linear polarization, m/. There 
is no clear relation between the two, and the obvious lack of 
correlation between linear and circular polarization is similar 
to the 5 GHz VLBA results of Homan et al. (2001), and tiie in- 
tegrated 5 GHz measurements from the ATCA by Rayner et al. 
(2000). The persistence of no trend at higher frequency sug- 
gests that the complicating factors (lack of resolution/Faraday 
effects) are to some degree self-similar and scale-up with ob- 
serving frequency. 

Although there is no correlation between linear and circular 
polarization for jet cores. Figure 5 shows that for quasars and 
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BL Lac objects, linear polarization almost always exceeds cir- 
cular polarization. This does not seem to be the case for radio 
galaxies, which generally have weak linear polarization. Brun- 
thaler et al. (200 1 ) have noted that the nearby very low luminos- 
ity radio sources Sgr A* and M8 1 * are circularly polarized but 
show no linear polarization. These authors have suggested that 
the ntc/mi ratio is distinctly different for low luminosity AGN 
{nic/mi » 1) and high luminosity AGN {rric/mi « 1). This 
relation was tested further by Bower et al. (2002), who failed to 
find additional LLAGN with circular polarization; however, if 
we consider the lowest redshift (and hence lowest luminosity) 
sources in our sample (which are admittedly quite different in 
character from the Sgr A* and M8 1 *), we find that nic / nii » I 
for M87 (1228+126) and 3C 84 (0316+413), the only two radio 
galaxies for which we can do the comparison. Both M87 and 
3 C 84 are known to have strong external Faraday screens which 
can depolarize the linear polarization (Wardle 1971; Zavala & 
Taylor 2002), and in fact, all the radio galaxies in the MOJAVE 
sample have very low levels of core linear polarization, likely 
due to external depolarization. Higher redshift sources may 
have smaller nic/mi ratios not because a different production 
mechanism for CP, but rather because either (1) they have less 
depolarization, and/or (2) we study them with less effective lin- 
ear resolution and therefore polarization from more of the jet 
contributes to what we see as the "core". If (2) is the case, it 
also naturally explains the lack of correlation between rric and 
mi in Figure 5. 

3.2. Circular Polarization of Jet Components 

Table 3 lists circular polarization measurements for all strong 
jet components in our sample with peak Stokes-I > 300 mJy 
beam"^. Of the 33 jet components that meet this criteria, we 
detect circular polarization in six (18%) of them at the 3ct 
level or higher significance. Five of these six detections are 
on the relatively nearby sources 3C84 (0316+413) and 3C273 
(1226+023) which are discussed in more detail in section §3.3, 
and we note that 3C 84 is a special case as the distinction be- 
tween the VLBI core and jet is not well defined and the cir- 
cular polarization measured is unusually strong (> 1%). For 
these reasons, we must be cautious about drawing conclusions 
about the circular polarization characteristics of jet components 
in general; however, it is clear that strong jet components can 
be circularly polarized, and the levels of circular polarization 
found in strong jet components are not clearly different from 
jet cores. 

3.3. Discussion of Selected Sources 

Here we discuss a few sources in detail, and we also com- 
pare our first epoch results with previous VLBA observations 
of circular polarization (HW99; Homan et al. (2001); Homan & 
Wardle (2003); Homan & Wardle (2004)) and from historical 
integrated measurements (1970s through early 1980s: Weiler 
& de Pater (1983), hereafter WdP83; Komesaroff et al. (1984), 
hereafter K84; AUer et al. (2003)). Because the latter are at 
lower frequency (typically 5 GHz and below) and are integrated 
measurements, we only compare to the historical circular polar- 
ization results in those cases where the source was repeatedly 
(three times or more) detected in circular polarization. 

Because there is a large overlap of our sample with the 
sources observed at 5 GHz with the VLBA on 1996 December 
15 by Homan et al. (2001), we Ust results from the latter study 
alongside our 15 GHz results in Table 5. We discuss some of 



the sources common to both of these studies in more detail be- 
low. The large number of limits from both programs makes the 
results in Table 5 difficult to interpret; however, it appears that 
CP of the core can vary with time and/or frequency. 

0316+413 (3C 84): The left-hand panels of Figure 6 show 
the central 5 milli-arcseconds from our March 2003 MOJAVE 
observations of 3C 84. The image reveals four distinct spots of 
circular polarization within this central region of the source. 
Without a clear component to identify as the core, we have 
taken spot "D" as the location of the VLBI core. Homan & War- 
dle (2004) found that even the locations corresponding to spots 
"C" and "B" were optically thick in December of 1997, so the 
precise definition of the jet core in this source is not clear With 
only one frequency, we don't know where the r ~ 1 surface is 
in 2003, so it is possible that all three spots of negative circular 
polarization: B, C, and D, may be part of an extended optically 
thick region; however, for the purposes of this paper, we have 
chosen to treat spots B and C as locations in the jet, along with 
spot A. 

Both spots D and C are similarly circularly polarized at the 
-0.4% level; however, the strength of the circular polarization 
increases sharply in spot B to -1.4%. A little further out in 
the jet, but at quite a different position angle, spot A is more 
than 2% circularly polarized with positive circular polarization. 
Homan & Wardle (2004) attribute this change in sign to the 
stark change in opacity they observe between spots B and A. 

The circular polarization we observe is very similar to that 
seen by Homan & Wardle (2004) in observations more than 

5 years earlier; however, spot D is entirely new and was not 
present in December of 1997. A version of the 15 GHz im- 
age from Homan & Wardle (2004) (convolved with a matching 
beam to our MOJAVE observations) appears in the right-hand 
panels of Figure 6. The 1997 image is registered with the 2003 
image by matching contours on the northern-most end of the 
central region shown here. The red dots in the right-hand pan- 
els of Figure 6 have the same relative positions as the blue stars 
which mark the locations of spots A, B, and C in 2003, but they 
have been shifted by 0.2 miUi-arcseconds to the North and by 
0.2 milli-arcseconds to the West. It is clear that the relative po- 
sitions of these circularly polarized spots has not changed sig- 
nificantly in 5+ years; however, they do appear to move outward 
together at a sub-luminal rate of roughly 0.06c, almost directly 
South-East. This speed is, of course, dependent on our align- 
ment of the images; however, it is similar to the speed found by 
Dhawan et al. (1998) at 43 GHz for two of the three component 
they studied within this region, although those two components 
were moving almost due South. 

Homan & Wardle (2004) found spot A to have = +3.2 ± 
0.1% which is somewhat larger than the nic = +2.3 ± 0.3% we 
measure five years later The intensity at the location of spot 
A has faded by a factor of three, from 1 .2 Jy beam"^ to 0.4 Jy 
beam~\ so presumably the region has become more optically 
thin. Multi-frequency observations are necessary to test if this 
is indeed the case, but if so, it is remarkable that this location 
has remained as strongly circularly polarized as it has. Homan 

6 Wardle (2004) argued that Faraday conversion is responsible 
for the high levels of circular polarization detected here. Fara- 
day conversion is strongly dependent on optical depth, and we 
would eventually expect a significant decline in circular polar- 
ization from this mechanism if this region of the jet is becom- 
ing increasingly optically thin (e.g., Wardle & Homan 2003). It 
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will be particularly interesting to follow the changes in spot A 
in future MOJAVE epochs. 

Homan & Wardle (2004) found spots B and C to have Wc = 
-0.7 ±0. 1% and rric = -0.6±0.1% respectively at 15 GHz. Spot 
B has essentially doubled in circular polarization percentage in 
5 years, and Spot C has not changed significantly within our 
uncertainties. Neither location has changed much in Stokes-/ 
intensity. 

If the alignment between the epochs is correct, the lack of 
large changes in the circular polarization distribution between 
spots A, B, and C is surprising given their motion. The mag- 
netic field order and particles responsible for producing this 
circular polarization must be in motion, however, this motion 
seems to have preserved the particular field characteristics re- 
sponsible for the strong levels of circular polarization we ob- 
serve. 

0430+052 (3C 120): We find a limit of |otc| < 0.29% in the 
core at 15 GHz in February of 2003. HW99 also found limits 
between 0.2 and 0.3% in the core at 15 GHz throughout 1996. 
Both signs of circular polarization (with a clear preference for 
positive CP) were detected intermittently in the historical inte- 
grated measurements (WdP83 and K84) at 8.9 GHz and below. 

0528+134: We find a limit of \ma\ < 0.18% in the core at 15 
GHz in February of 2003. HW99 detected multiple epochs of 
significant circular polarization at about the ntc = -0.5% level in 
the core at 15 GHz during 1996. During 1996, 0528+134 was 
undergoing a large outburst, peaking with an integrated core 
flux 4.6 times larger than observed in February 2003. 

0607-157: We find a hmit of {m,] < 0.21% in the core at 
15 GHz in March of 2003. Homan & Wardle (2003) measured 
only Me = -0.23 ± 0.12% in the core at 15 GHz in January of 
1998, but they simultaneously detected nic = -0.65 ±0.07% in 
the core at 8 GHz. Homan et al. (2001) detected circular po- 
larization in core of nic = -0.18 ±0.05 % at 5 GHz thirteen 
months earlier in December of 1996. Aller et al. (2003) found 
strong circular polarization, iric = -0.23 ± 0.04% at 4.8 GHz in 
integrated measurements with the UMRAO during the period 
2001-2002. 

0735+178: We measure nic = -0.30 ± 0. 1 1 % in the core at 
15 GHz in November of 2002. HW99 did not detect core cir- 
cular polarization at 15 GHz at any epoch during 1996 with a 
typical hmit of \mc\ < 0.4%. 

0851+202 (OJ287): We measure = -0.20 ± 0.08 % in the 
core at 15 GHz in October of 2002. HW99 did not detect core 
circular polarization at 15 GHz at any epoch during 1996, with 
a typical limit of \mc\ < 0.2%. Homan et al. (2001) also found a 
limit of |mc I < 0. 10% for the core circular polarization at 5 GHz 
in December of 1996. In older integrated measurements during 
the period 1978-1983, Aller et al. (2003) found strong CP at 4.8 
and 8.0 GHz with nic = -0.37 ± 0. 12% and nic = -0.22 ± 0.03% 
respectively. 

0923+392 (4C 39.25): We find a limit of |m,| < 0.18% in the 
core at 15 GHz in October of 2002. In older integrated measure- 
ments during flie period 1978-1983, Aller et al. (2003) found 
strong CP at 4.8 GHz with nic = -0.21 ± 0.05%. 
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1127-145: We measure nic = -0.24 ± 0.10 % in the core at 
15 GHz in October of 2002. Negative circular polarization was 
repeatedly detected at 5 GHz and below in the historical inte- 
grated measurements cataloged by WdP83 and K84. 



1222+216: We find a hmit of \mc\ < 0.35% in flie core at 15 
GHz in February of 2003. HW99 also found limits between 0.1 
and 0.2% in the core at 15 GHz throughout 1996. 



1226+023 (3C273): Figure 7 shows flie first 10 miUi- 
arcseconds of the jet of 3C 273 at 15 GHz in October of 2002. 
The top panel is Stokes I, the middle panel is hnear polariza- 
tion, and the bottom panel is circular polarization. Positions 
where circular polarization measurements or limits were taken 
are marked and numbered on the Stokes 1 map. Figure 8 shows 
the fractional circular polarization at these locations as a func- 
tion of distance from the jet core. The absolute fractional level 
of the circular polarization appears to decrease as distance from 
the core increases. This decrease could be tied to decreasing 
opacity along the jet or a decrease in poloidal magnetic field 
strength. 

The core of 3C 273 was found by HW99 to have nic - -0.5% 
at three epochs during 1996. The circular polarization appeared 
with the emergence of a new jet component. The core was very 
opaque at this time, and HW99 showed that the appearance of 
the CP was tied to the flattening of the core's spectral index 
due to the emergence of the new component. As cataloged in 
WdP83 and K84, 3C 273 was repeatedly detected in integrated 
measurements during the 1970s and early 80s at 8 GHz and 
below where it had a consistently negative sign of circular po- 
larization. 



1228+126 (M87): Figure 9 shows our 15 GHz VLBA image 
of the jet of M87 in February of 2003. No Unear polarization is 
detected, but strong circular polarization, ntc = —0.49 ± 0.10% 
is seen at the peak of the core. 

M87 was studied at the VLA in AprU of 2000 by Bower et al. 
(2002), who found no significant circular polarization at 8 GHz, 
me = 0.01±0.10; however, they did find linear polarization at 
the 1.7% level. These results are the opposite of our finding 
of strong circular polarization and no hnear polarization. It is 
important to consider that their VLA results include contribu- 
tions from the jet weU beyond the region imaged by the VLBA, 
and the linear polarization they observe must come from be- 
yond the depolarizing screen which covers most of the VLBA 
scale jet (Zavala & Taylor 2002). The fractional circular po- 
larization they measure may be diluted by emission beyond the 
VLBA scale, or the circular polarization may simply vary with 
frequency or time. Regardless, the differing linear polariza- 
tions observed between the VLBA and VLA scale images un- 
derscores our point from §3.1 that at least some of the lack of 
correlation between Unear and circular polarization for jet cores 
could be due to the two quantities being produced in different 
regions of the unresolved jet core. 

M87 has been detected intermittently in integrated circular 
polarization (with both signs) at 5 GHz and below in the histor- 
ical observations cataloged by WdP83 and K84. 
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1253-055 (3C279): Figure 10 shows the first 5 milh- 
arcseconds of the jet of 3C279 at 15 GHz in November of 
2002. The top panel is Stokes I, the middle panel is Unear po- 
larization, and the bottom panel is circular polarization. Posi- 
tions where circular polarization measurements or limits were 
taken are marked and numbered on the Stokes I map. We detect 
strong circular polarization of rric - +0.30 ± 0.08% only on the 
jet core. The sign and magnitude of the circular polarization 
is consistent with five epochs of CP detections by HW99 at 15 
GHz during 1996. At 5 GHz, Homan et al. (2001) did not de- 
tect circular polarization with a limit of \mc\ < 0.11%. Histor- 
ically, 3C 279 has repeatedly and consistently shown a positive 
sign for its circular polarization in integrated measurements at 
8 GHz and below during the 1970s and early 80s (WdP83 and 
K84); however, Aller et al. (2003) found nic = -0.17 ±0.01% 
at 4.8 GHz with the UMRAO during 2001-2002. This change 
in sign at the lower frequency is particularly interesting as no 
change has yet occurred at 15 GHz. 

1308-1-326: We find a limit of \mc | < 0.21 % in the core at 15 

GHz in November of 2003. HW99 also found limits between 
0.1 and 0.2% in the core at 15 GHz throughout 1996. 

1510-089: We measure = +0.20 ± 0.09 % in the core at 15 
GHz in November of 2002. HW99 did not detect circular po- 
larization in the core at 15 GHz during 1996: three epochs had 
limits of 0.1 to 0.2%, and for two epochs they reported peak 
CP values of +0.2 ± 0.2 %. Homan et al. (2001) find a limit of 
|otc| < 0.24% in the core at 5 GHz in December of 1996. 

1641+399 (3C345): We find a limit of [md < 0.19% in 
the core at 15 GHz in November of 2002. Homan & War- 
die (2003) detected circular polarization in the core at 15 GHz, 
iric =-0.38 ±012%, in January of 1998. The same location was 
optically thick at 8 GHz where they did not detect any circular 
polarization in simultaneous measurements. 3C 345 was de- 
tected intermittently (and with both signs) at 8 GHz and below 
in the historical integrated measurement compiled by WdP83. 

1730-130: We find a Umit of \mc\ < 0.17% in the core at 
15 GHz in October of 2002. The historical integrated measure- 
ments (WdP83 and K84) show several intermittent detections 
of positive circular polarization at 5 GHz and below. 

1749+096: We find a limit of \mc \ < 0. 14% in the core at 15 
GHz in May of 2002. HW99 also found typical Umits of 0.2% 
in the core at 15 GHz throughout 1996. 

1928+738: We measure = +0. 19 ± 0.09 % in the core at 
15 GHz in June of 2002. HW99 also measured circular polar- 
ization in the core 15 GHz at the limit of their sensitivity. They 
found nic = +0.3 ±0.1 % in two epochs, nic = +0.2 ± 0.1 % in 
one epoch, and limits of 0. 1 and 0.2% in their other two epochs 
during 1996. 

2134+004: Figure 11 shows our observations of 2134+004 
at 15 GHz in May of 2002. The top panel is Stokes 1, the 
middle panel is linear polarization, and the bottom panel is cir- 
cular polarization. Positions where circular polarization mea- 
surements or limits were taken are marked and numbered on 



the Stokes 1 map. We detect strong circular polarization of 
nic = +0.35 ±0.12% on the jet core, we measure a similar value 
at location number 1 of me = +0.29 ±0. 12%, and we find a Umit 
on the CP at location 2 to be jm^ | < 0.34%. 

In recent integrated measurements from 2001-2002, Aller et 
al. (2003) found strong circular polarization at 4.8 GHz with 
nic = +0.62 ±0.06%. 2134+004 was also detected repeatedly 
in the historical integrated measurements at 5 GHz and below, 
always with positive circular polarization (WdP83 and K84). 



2145+067: We find a Umit of Im^ | < 0.26% in the core at 15 
GHz in March of 2003. In recent integrated measurement from 
2001-2002, AUer et al. (2003) measured nic = -0.50 ± 0.07% at 
4.8 GHz with the UMRAO. The historical integrated measure- 
ments (WdP83 and K84) show consistent, repeated detections 
of negative circular polarization at 5 GHz and below. Given this 
historical record of detection, it is interesting that the 'limit' re- 
ported here came from measurements that had a significance of 
1 .98(7, just barely less than the 2a cutoff we had for defining 
measurements versus limits. If we allow 1.98(t as a measure- 
ment, then our fractional level for the core is nic = -0.17 ± 0.09 
%. 



2200+420: We find a Umit of \mc\ < 0.21% in the core at 
15 GHz in June of 2002. WdP83 reported intermittent detec- 
tions of both signs of circular polarization at frequencies below 
5 GHz in the historical integrated measurements. 



2230+114 (CTA102): We find a Umit of |mc | < 0. 1 9% in the 
core at 15 GHz in October of 2002. Homan et al. (2001) mea- 
sured nic = -0.09 ± 0.04% in the core at 5 GHz in December of 
1996. CTA102 was intermittently detected (with both signs of 
CP) in the historical integrated measurements (WdP83; K84) at 
5 GHz and below. 



2251+158 (3C 454.3): Figure 12 shows the first 4 milU- 
arcseconds of the jet of 3C 454.3 at 15 GHz in March of 
2003. Positions where circular polarization measurements or 
limits were taken are marked and numbered on the Stokes I 
map. The first jet component (position 1) appears sUghtly 
more strongly circularly polarized, = +0.38 ±0.10%, than 
the core, nic = +0.23 ± 0.10%; however, an alternative calibra- 
tion assuming no circular polarization for calibrating the phases 
(see §2.2) spreads the circular polarization between the core and 
component 1 more evenly with both having ntc — 0.3%. 

Homan et al. (2001) measured nic = +0.08 ±0.04% in the 
core of 3C 454.3 at 5 GHz in December of 1996, and the source 
was repeatedly detected with positive circular polarization at 8 
GHz and below in the integrated historical measurements com- 
piled byWdP83andK84. 



2345-167: We detect nic = +0.40 ± 0.11% in the core at 15 
GHz in October of 2002. Homan et al. (2001) found a limit 
of I mc I < 0.11% in the core at 5 GHz in December of 1996. 
WdP83 and K84 reported intermittent detections (but always 
with positive sign) of circular polarization in integrated mea- 
surements at frequencies of 5 GHz and below. 
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4. CONCLUSIONS 

In this paper, we report first-epoch circular polarization re- 
sults for each of the 133 AGN in the MOJAVE program to 
monitor the VLBA structure and polarization of a flux limited 
sample of AGN jets at 15 GHz. We found strong circular po- 
larization (> 0.3%) in approximately 15% of our sample. The 
circular polarization we detected was usually associated withjet 
cores, typically in the range of 0.3-0.5% of the local Stokes- 
/. Our ability to detect circular polarization on jet components 
was limited by their relatively weak flux density; however, we 
did find a few strong jet components to be circularly polarized. 

We found no significant correlations between fractional cir- 
cular polarization of jet cores and optical source type, redshift, 
EGRET detections, linear polarization, or a number of other pc- 
scale jet properties listed in Table 4. While there was no clear 
relation between linear and circular polarization of jet cores, 
quasars and BL Lac objects consistently had more linear than 
circular polarization, and the two galaxies with detected circu- 
lar polarization had much more circular than hnear polarization. 
We suggested that higher redshift sources may have smaller 
mc/mi ratios, not because they produce circular polarization 
differently from low redshift sources, but rather because either 
(1) they have less depolarization, and/or (2) we study them with 
less effective linear resolution and therefore polarization from 
more of the jet contributes to what we see as the "core". 

The lack of clear correlations between circular polarization 
and other source properties is intriguing, particularly given 
that some sources, 3C 273 and 3C 279 in particular, show con- 
sistently high levels of circular polarization over years and 
decades of observation, even though they are both known to 
be highly variable with regular (~ annual) ejections of new jet 
components. It seems highly unhkely that these sources would 
be circularly polarized in a very similar fashion from epoch to 



epoch due only to stochastic factors (such as a chance favorable 
magnetic field configuration). It may be that we simply detect 
too few sources in circular polarization to find subtle trends be- 
tween circular polarization and other source properties, or per- 
haps, even with the VLBA we are not resolving the core region 
weU enough to see where the circular polarization is being gen- 
erated. 

In addition to looking for overall trends in the circular polar- 
ization of our sample, we also examined the structure of a few 
sources in detail. 3C 84 was particularly interesting, showing 
the strongest circular polarization in our sample, up to 2% in 
one location. Our image revealed essentially the same structure 
as Homan & Wardle (2004) found in data taken five years ear- 
lier; however, a new circularly polarized spot appeared in our 
data, and the older spots of circular polarization, while main- 
taining the same relative positions, appear to have moved out- 
ward in the jet at a rate of approximately 0.06c. 3C 273 also 
showed multiple locations of circular polarization in its jet, and 
we were able to plot circular polarization measurements for sev- 
eral strong jet components. This plot revealed that, in 3C 273, 
fractional circular polarization gets weaker at greater distances 
from the jet core. Multi-frequency observations will be neces- 
sary for both sources so we can use opacity measurements to 
help understand the trends of circular polarization with jet po- 
sition. 

The authors wish to acknowledge the contributions of the 
other members of the MOJAVE team: Hugh and Margo 
Aller, Tigran Arshakian, Marshall Cohen, Matthias Kadler, 
Ken Kellermann, Yuri Kovalev, Andrei Lobanov, Eduardo Ros, 
Rene Vermeulen, and Tony Zensus. This work has been sup- 
ported by Denison University, the National Radio Astronomy 
Observatory, and NSF grant 0406923-AST. This research was 
also supported by an award from Research Corporation. 
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Circular Polarization at Stokes / Peak. 
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First-Epoch Circular Polarization in MOJAVE 
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PHL 5225 


1.562 


Q 


2003 Mar 29 


1.88 


0.34 


+0.67 ± 1.71 


< 0.18 




2230+114 


CTA 102 


1.037 


Q 


2002 Oct 09 


2.32 


1.67 


+0.12 ±2.22 


< 0.19 




2243-123 




0.63 


Q 


2002 Oct 09 


1.03 


1.32 


+0.97 ± 1.27 


< 0.25 




2251 + 158 


3C 454.3 


0.859 


Q 


2003 Mar 29 


3.25 


1.05 


+7.35 ±3.32 


+0.23 ±0.10 


2.2 


2331+073 






U 


2002 Jun 15 


0.96 


3.38 


+1.26 ±1.06 


<0.24 




2345-167 




0.576 


Q 


2002 Oct 09 


1.39 


0.77 


+5.55 ± 1.48 


+0.40 ±0.11 


3.8 


2351+456 


4C +45.51 


1.986 


Q 


2002 May 31 


1.39 


2.09 


+2.83 ±1.45 


<0.31 





Note. — Columns are as follows: (1) lAU Name (B1950); (2) Other Name; (3) Redshift (see Lister & Homan 2005, for references); (4) 
Optical classification where Q = quasar, B = BL Lac object, G = active galaxy, and U = unidentified; (5) Epoch of first observation in MOJAVE 
program; (6) Stokes / at the map peak; (7) Fractional linear polarization of the Stokes / peak; (8) Maximimi circular polarization measured 
within two pixels of the Stokes / peak; (9) Fractional circular polarization corresponding to Stokes-V in column (8); (10) Significance of the 
circular polarization measurement: | | / Am^ . 
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Table 2 



Circular Polarization on Off-Peak Jet Cores. 



lAV 


R 


e 


leore 


mi 


Vco„ 


trie 




Name 


(marcsec) 


(deg.) 


(Jybeam"') 


(%) 


(mjy beam"') 


(%) 


a 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


0212 + 735 


0.67 


-63.4 


0.79 


5.97 


+0.21 ±0.90 


<0.23 




0316+413 


0.28 


45.0 


2.26 


<0.05 


-8.61 ±3.15 


-0.38 ±0.14 


2.7 


0738 + 313 


3.71 


-4.6 


0.33 


< 0.32 


+0.37 ± 1.12 


<0.68 




0923 + 392 


2.75 


-79.5 


0.19 


< 0.78 


-0.95 ±0.81 


< 0.92 




1213-172 


0.50 


-53.1 


0.53 


1.84 


+0.86 ±0.71 


<0.30 




1226 + 023 


0.72 


56.3 


2.78 


0.62 


-18.77 ±2.71 


-0.67 ±0.10 


6.9 


1548 + 056 


0.63 


161.6 


1.25 


1.51 


-0.20 ±1.29 


<0.21 




2128-123 


1.12 


26.6 


0.43 


<0.20 


-0.20 ±0.51 


<0.23 





Note. — Columns are as follows: (1) lAU Name (B1950); (2) Radial position of the core relative to the map 
peak; (3) Structural position angle of the core relative to the map peak; (4) Peak Stokes / at the core position; (5) 
Fractional linear polarization at the core position; (6) Maximum circixlar polarization measured within two pixels 
of the core; (7) Fractional circular polarization corresponding to Stokes V in coluinn(6); (8) Significance of the 
circular polarization measiffement: |mc I / Am^. 



Table 3 

Circular Polarization on Strong Jet Components. 



lAU 


R 


e 




mi 








Name 


(marcsec) 


(deg.) 


(Jybeam"') 


(%) 


(mjybeam"') 


(%) 


cr 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


0212 + 735 


0.67 


116.6 


1.25 


10.90 


+2.46 ±1.35 


<0.31 




0224 + 671 


0.41 


14.0 


0.60 


1.17 


+0.12 ±0.78 


<0.26 




0316+413 


1.00 


-143.1 


2.00 


<0.06 


-8.43 ±2.82 


-0.42 ±0.14 


3.0 


0316+413 


1.77 


-137.3 


0.65 


< 0.17 


-8.93 ± 1.13 


-1.37±0.17 


7.9 


0316 + 413 


1.94 


-168.1 


0.34 


< 0.34 


+7.56 ±0.84 


+2.25 ±0.25 


9.0 


0552 + 398 


0.85 


-69.4 


0.33 


<0.37 


-0.34 ±0.48 


<0.29 




0738 + 313 


3.71 


175.4 


0.51 


0.42 


-1.79 ±1.65 


<0.68 




0754+100 


1.22 


9.5 


0.31 


4.03 


-0.13 ±0.39 


< 0.25 




0851+202 


1.25 


-118.6 


0.37 


7.73 


+0.61 ±0.52 


<0.30 




0923 + 392 


2.75 


100.5 


4.33 


2.05 


+4.09 ±3.76 


<0.18 




0945 +408 


0.94 


122.0 


0.33 


i.n 


+0.75 ±0.44 


<0.36 




1038 + 064 


0.71 


171.9 


0.49 


1.29 


-0.81 ±0.56 


<0.28 




1127-145 


0.51 


101.3 


0.57 


7.30 


-0.35 ±0.63 


< 0.22 




1213-172 


0.50 


126.9 


1.11 


4.18 


+0.61 ± 1.39 


< 0.25 




1226 + 023 


0.72 


-123.7 


8.99 


3.74 


-40.28 ± 8.07 


-0.45 ± 0.09 


5.0 


1226 + 023 


2.56 


-120.6 


0.46 


6.33 


-2.59 ±1.18 


-0.56 ±0.26 


2.2 


1226 + 023 


3.55 


-130.4 


1.28 


7.08 


-4.88 ±1.59 


-0.38 ±0.12 


3.1 


1226 + 023 


4.53 


-131.4 


0.82 


4.20 


-1.00 ± 1.33 


< 0.32 




1226 + 023 


9.75 


-115.5 


1.49 


9.97 


+0.10± 1.73 


< 0.23 




1253 - 055 


1.36 


-126.0 


1.04 


24.91 


-0. 10 ± 1.65 


< 0.32 




1253 - 055 


2.83 


-122.0 


0.52 


38.68 


+0.38 ± 1.47 


<0.56 




1253 - 055 


4.99 


-122.7 


0.54 


10.93 


-1.40 ± 1.47 


< 0.55 




1548 + 056 


0.63 


-18.4 


1.66 


3.24 


-0.38 ± 1.68 


< 0.20 




1641 + 399 


1.10 


-90.0 


1.53 


4.14 


+1.05 ± 1.66 


< 0.22 




1641+399 


2.31 


-95.0 


0.33 


3.43 


-0.77 ± 0.96 


< 0.58 




1928 + 738 


0.98 


156.0 


0.52 


1.58 


+1.09 ±0.61 


< 0.33 




2128-123 


1.12 


-153.4 


1.78 


1.72 


+3.86 ±1.63 


+0.22 ± 0.09 




2131-021 


0.76 


113.2 


0.61 


1.92 


+0.46 ±0.63 


<0.21 




2134 + 004 


0.80 


-90.0 


0.77 


8.54 


+2.21 ±0.94 


+0.29 ±0.12 




2134 + 004 


2.02 


-98.5 


1.33 


4.99 


+2.99 ±1.56 


< 0.34 




2145 + 067 


0.71 


135.0 


1.83 


3.01 


+0.59 ± 1.76 


< 0.19 




2223-052 


0.85 


45.0 


0.63 


6.03 


+1.20 ±0.76 


<0.31 




2230 + 114 


1.20 


138.4 


0.32 


7.05 


+0.11 ±0.45 


<0.28 




2243-123 


1.60 


-3.6 


0.43 


2.30 


+0.60 ±0.57 


<0.27 




2251 + 158 


0.80 


-90.0 


3.11 


1.10 


+11.94±3.20 


+0.38 ±0.10 


3.7 



Note. — Columns are as follows: (1) lAU Name (B 1950); (2) Radial position of the component relative to the jet 
core; (3) Structural position angle of the component relative to the jet core; (4) Peak Stokes / of the component; (5) 
Fractional linear polarization of the component; (6) Maximtmi circular polarization measured within two pixels of 
the component peak; (7) Fractional circular polarization corresponding to Stokes-V in column (6); (8) Significance 
of the circular polarization measurement: \mc\/Amc. 
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Table 4 

Kendall-Tau Test for Correlations with Fractional Core CP. 



Property p 



Apparent Optical V Magnitude 0.07 

Optical V Magnitude 0.78 

Redshift 0.33 

Galactic Latitude, 0.28 

Luminosity Distance 0.33 

Angular Size Distance 0.10 

Maximum Jet Speed, (3app 0.68 

Total 15 GHz Flux 0.57 

Maximum Total 15 GHz Flux (any epoch) 0.44 

Total 15 GHz Luminosity 0.38 
Maximum Total 15 GHz Luminosity (any epoch) 0.66 

15 GHz Flux of Core 0.94 

15 GHz Luminosity of Core 0.30 

15 GHz Flux of Jet 0.06 

15 GHz Luminosity of Jet 0.92 

Core to Jet Flux Ratio 0.05 

Core to Total Flux Ratio 0.04 

Total Fractional Linear Polarization 0.44 

Fractional Linear Polarization of Jet 0.46 

Fractional Linear Polarization of Core 0.19 



Note. — The /) value indicates the probability (1 = 100%) 
that no correlation exists between the (absolute) fractional 
circular polarization of the core, \mc\, and the listed prop- 
erty. References for optical data are given in Lister & Homan 
(2005). Proper motion data is from Kellermann et al. (2004). 
Radio data is from Lister & Homan (2005). 



Table 5 

Comparison of Core CP at 5 GHz and 15 GHz. 



IAU Name 


5 GHz rric (%) 


15 GHz m. 


(%) 


0215+015 


< 0.21 


+0.28 ± 0. 


12 


0336-019 


+0.19 ±0.06 


<0.18 




0403-132 


<0.16 


<0.18 




0420-014 


<0.11 


<0.16 




0605-085 


<0.10 


< 0.26 




0607-157 


-0.18 ±0.05 


<0.21 




0851+202 


<0.10 


-0.20 ±0. 


08 


0906+015 


<0.43 


<0.19 




1055+018 


<0.13 


+0.32 ±0. 


09 


1156+295 


< 0.11 


-0.27 ± 0. 


09 


1253-055 


< 0.11 


+0.30 ± 0. 


08 


1334-127 


-0.12 ±0.05 


+0.29 ±0. 


10 


1413+135 


<0.18 


<0.34 




1502+106 


<0.15 


<0.20 




1510-089 


< 0.24 


+0.20 ± 0. 


09 


1546+027 


-0.25 ± 0.05 


< 0.24 




2201+171 


+0.24 ±0.06 


+0.19 ±0. 


09 


2223-052 


+0.18±0.04 


<0.22 




2227-088 


<0.12 


<0.18 




2230+114 


-0.09 ± 0.04 


< 0.19 




2243-123 


-0.14 ±0.04 


< 0.25 




2251+158 


+0.08 ±0.04 


+0.23 ±0. 


10 


2345-167 


<0.11 


+O.40±0. 


11 



Note. — The 15 GHz data are from Tables 1 and 2 
of this paper. The 5 GHz data (epoch 1996 December 
15) are from Homan et al. (2001) 
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Fig. 1 . — Gain calibration uncertainty for circular polarization (as a percentage of Stokes-/) as a function of Stokes-/ intensity at the map peak for tiie 133 sources 
in our sample. 
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Fig. 2. — Circular polarization (color-scale given in mjy beam ) superimposed on Stokes-/ contours for 1213—172 (top panels) and 2243—123 (bottom panels). 
The / contours start at 1 and 2 mJy beam"' respectively for 1213—172 and 2243-123 and increase in steps of X\/2. The FWHM dimensions of the restoring beams 
are indicated by a cross in the lower left comer of each panel. Left hand panels (a,c) show the results of straight gain transfer calibration which has a combination of 
real circular polarization plus phase errors. Right hand panels (b,d) show the additional effect of a round of phase self-calibration assuming no circular polarization 
in the data. For both sources, the additional phase calibration has corrected a phase shift in the data between the left and right hand data, and reveal that the source 
has no real circular polarization. 
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Fig. 3. — Circulai' polaiization (color-scale given in mJy beam ) superimposed on Stokes-/ contours for 3C 84 (top panels) and 3C 273 (bottom panels). The / 
contours for both sources start at 10 mJy beam"' and increase in steps of X \/2. The FWHM dimensions of the restoring beams are indicated by a cross in the lower 
left comer of each panel. Left hand panels (a,c) show the results of straight gain transfer calibration which has a combination of real circular polarization plus phase 
errors. Right hand panels (b,d) show the additional effect of a round of phase self-calibration assuming no circular polarization in the data. For 3C 84, panel (b), the 
additional phase calibration has corrected the phase errors and revealed the circular polarization of the source in greater detail, although the main features appear 
quite similar to panel (a). For 3C 273, panel (d), the additional phase calibration has corrected the phase error and reveals that the real circular polarization in the 
source is on the jet core and not on the jet component. Note that the images shown here for 3C273 are from our second MOJAVE epoch on this source (August 28, 
2003), but they are included here as an illustrative example. 
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Fig. 4. — Histograms of the absolute value of the fractional circular polarization at the location of the VLB! jet core with each of the bottom three panels 
con'esponding to a different optical classification. The significance of each measurement is indicated by open, hash, or solid fill styles for < 2cr, 2-3(7, and > 3cr 
values respectively. Note that no limits appear on these histograms. For sources with limits in Table 1 or 2, we have substituted the fractional circular polarization 
at the core location calculated from the V and / values reported in the table. 
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Fig. 5. — Plot of the absolute value of the fractional circular polarization versus fractional linear polarization for jet cores. Black points represent quasars, red 
points are BL Lac objects, blue points are galaxies, and unidentified objects appear in green. Solid squares are of > 3a significance, open squares are of 2-3(7 
significance, and arrows represent upper limits. The dashed line indicates \mc\ = mi. 
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Fig. 6. — Circular polarization from the central 5 milli-arcseconds of 3C 84. The left-hand panels show our March 2003 observations, and the right-hand panels 
show the observations by Homan & Wardle (2004) from December of 1997. The images from both epochs have been convolved with an identical beam of 0.69x0.55 
mas X mas, plotted as a cross in the lower left hand corner of the plot. The top images are circular polarization contours beginning at ±2 mjy beam"' and increasing 
in steps of X \/2. Dashed contours represent negative circular polarization. A single Stokes-/ contour shows the registration of the CP within this region. The bottom 
images are the Stokes-/ contours (beginning at ±10 mJy beam"' and increasing in steps of X \f2) from the same region. The four circularly polarized components 
from March 2003 are labelled with letters and their positions are marked with red stars. The three blue dots on the right hand image have the same relative positions 
as the red stars corresponding to components A-C, only they have been shifted by 0.2 milli-arcseconds to the west and by 0.2 milli-arcseconds to the north. 
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Fig. 7. — Stokes-/ contours (top), linear polarization (middle), and circular polarization (bottom) images of the inner 10 milli-arcseconds of the jet of 3C273 in 
October of 2002. A single Stokes-/ contour appears around the linear and circular polarization images to show registration. The Stokes-/ and linear polarization 
contours begin at 20 mjy beam"' , and the circular polarization contours begin at ±2 mjy beam"' . All contours increase in steps of x \/2. In this image, negative 
circular polarization is indicated by solid contours. The FWHM dimensions of the common restoring beam is indicated by a cross in the lower left comer of the 
image. Electric vector position angles (EVPA) for the hnear polarization are indicated by tick marks in the middle image. The locations of circular polarization 
measurements are indicated on the Stokes-/ map. 
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Fig. 8 . — Plot of fractional circular polarization as a function of distance from the jet core for 3C 273. Note that the circular polarizaiton in 3C 273 predominantly 
has a negative sign and becomes weaker (going to smaller absolute values) at greater distances from the core. For the purposes of this plot, we have plotted all the 
values with their corresponding error bars, including those of < 2a significance. 
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Fig. 9. — Stokes-/ contours (top), linear polarization (middle), and circular polarization (bottom) images of the VLBI jet of M87 in February of 2003. A single 
Stokes-/ contour appears around the linear and circular polarization images to show registration. All contours begin at ±1 mjy beam"' and increase in steps of 
X \/2. Negative circular polarization is indicated by dashed contours. The FWHM dimensions of the common restoring beam is indicated by a cross in the lower 
left comer of the image. 
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Fig. 10. — Stokes-/ contours (top), linear polarization (middle), and circular polarization (bottom) images of tlie inner 5 miUi-arcseconds of the jet of 3C279 in 
November of 2002. A single Stokes-/ contour appears around the linear and circular polarization images to show registration. The Stokes-/ and linear polarization 
contours begin at 20 mjy beam"', and the circular polarization contours begin at ±5 mjy beam"'. All contours increase in steps of x^/l. Negative circular 
polarization is indicated by dashed contours. The FWHM dimensions of the common restoring beam is indicated by a cross in the lower left comer of the 
image. Electric vector position angles (EVPA) for the linear polarization are indicated by tick marks in the middle image. The locations of circular polarization 
measurements are indicated on the Stokes-/ map. 
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Fig. 11. — Stokes-7 contours (top), linear polarization (middle), and circular polarization (bottom) images of 2134+004 in May of 2002. A single Stolces-/ 
contour appears around the linear and circular polarization images to show registration. The Stokes-/ and linear polarization contours begin at 5 mjy beam~\ and 
the circular polarization contours begin at it 1 mJy beam"' . All contours increase in steps of X \/2. Negative circular polarization is indicated by dashed contours. 
The FWHM dimensions of the common restoring beam is indicated by a cross in the lower left comer of the image. Electric vector position angles (EVPA) for the 
linear polarization are indicated by tick marks in the middle image. The locations of circular polarization measurements are indicated on the Stokes-/ map. 
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Fig. 12. — Stokes-/ contours (top), linear polarization (middle), and circular polarization (bottom) images of the inner 4 milli-arcseconds of the jet of 3C 454.3 
in March of 2003. A single Stokes-/ contour appears around the linear and circular polarization images to show registration. The Stokes-/ and linear polarization 
contours begin at 5 mJy beam"' , and the circular polarization contours begin at ±3 mJy beam"' . All contours increase in steps of X Negative circular 
polarization is indicated by dashed contours. The FWHM dimensions of the common restoring beam is indicated by a cross in the lower left corner of the 
image. Electric vector position angles (EVPA) for the linear polarization are indicated by tick marks in the middle image. The locations of circular polarization 
measurements are indicated on the Stokes-/ map. 



